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Abstract. The current status of leptonic jet models for gamma-ray blazars is reviewed. 
Differences between the quasar and BL-Lac subclasses of blazars may be understood 
in terms of the dominance of different radiation mechanisms in the gamma-ray regime. 
Spectral variability patterns of different blazar subclasses appear to be significantly 
different and require different intrinsic mechanisms causing gamma-ray flares. As ex- 
amples, recent results of long-term multiwavelength monitoring of PKS 0528-f 134 and 
Mrk 501 are interpreted in the framework of leptonic jet models. A simple quasi- 
analytic toy model for broadband spectral variability of blazars is presented. 



INTRODUCTION 

Recent high-energy detections and simultaneous broadband observations of 
blazars, determining their spectra and spectral variability, are posing strong con- 
straints on currently popular jet models of blazars. 66 blazars have been detected 
by EGRET at energies above 100 MeV [1], the two nearby high-frequency peaked 
BL Lac objects (HBLs) Mrk 421 and Mrk 501 are now multiply confirmed sources 
of multi-GeV - TeV radiation [2-5], and the TeV detections of PKS 2155-314 [6] 
and lES 2344-1-514 [7] are awaiting confirmation. Most EGRET-detected blazars 
exhibit rapid variability [8], in some cases on intraday and even sub- hour (e. g., [9]) 
timescales, where generally the most rapid variations are observed at the highest 
photon frequencies. 

The broadband spectra of blazars consist of at least two clearly distinct spec- 
tral components. The first one extends in the case of flat-spectrum radio quasars 
(FSRQs) from radio to optical/UV frequencies, in the case of HBLs up to soft and 
even hard X-rays, and is consistent with non-thermal synchrotron radiation from 
ultrarelativistic electrons. The second spectral component emerges at 7-ray ener- 



gies and peaks at several MeV - a few GeV in most quasars, while in the case of 
some HBLs the peak of this component appears to be located at TeV energies. 

The bolometric luminosity of EGRET-detected quasars and some low-frequency 
peaked BL Lac objects (LBLs) during flares is dominated by the 7-ray emission. 
If this emission were isotropic, it would correspond to enormous luminosities (up 
to ~ 10^^ erg s~^) which, in combination with the short observed timescales (im- 
plying a small size of the emission region) would lead to a strong modiflcation of 
the emissivity spectra by 77 absorption, in contradiction to the observed smooth 
power-laws at EGRET energies. This has motivated the concept of relativistic 
beaming of radiation emitted by ultrarelativistic particles moving at relativistic 
bulk speed along a jet (for a review of these arguments, see [10]). While it is 
generally accepted that blazar emission originates in relativistic jets, the radiation 
mechanisms responsible for the observed 7 radiation are still under debate. It is 
not clear yet whether in these jets protons are the primarily accelerated particles, 
which then produce the 7 radiation via photo-pair and photo-pion production, fol- 
lowed by 7r° decay and synchrotron emission by secondary particles (e. g., [11]), or 
electrons (and positrons) are accelerated directly and produce 7-rays in Compton 
scattering interactions with the various target photon fields in the jet [12-15]. 

In this review, I will describe the current status of blazar models based on lep- 
tons (electrons and/or pairs; in the following, the term "electrons" refers to both 
electrons and positrons) as the primary constituents of the jet which are responsible 
for the 7-ray emission. Hadronic jet models are discussed in a separate paper by J. 
Rachen [16]. In Section 2, I will give a description of the model and discuss the dif- 
ferent 7-ray production mechanisms and their relevance for different blazar classes. 
In Section 3, I will review recent progress in understanding intrinsic differences be- 
tween different blazar classes and present state-of-the-art model calculations, using 
a leptonic jet model, to undermine the general theoretical concept. In Section 4, 
I will discuss broadband spectral variability of individual blazars and their inter- 
pretation in the framework of leptonic jet models. A simple quasi-analytical toy 
model for blazar broadband spectral variability will be presented in section 5. 

MODEL DESCRIPTION AND RADIATION 
MECHANISMS 

The basic geometry of leptonic blazar jet models is illustrated in Fig. 1. At the 
center of the AGN, an accretion disk around a supermassive, probably rotating, 
black hole is powering a relativistic jet. Along this pre-existing jet structure, occa- 
sionally blobs of ultrarelativistic electrons are ejected at relativistic bulk velocity. 

The electrons are emitting synchrotron radiation, which will be observable at 
IR - UV or even X-ray frequencies, and hard X-rays and 7-rays via Compton 
scattering processes. Possible target photon fields for Compton scattering are the 
synchrotron photons produced within the jet (the SSC process, [12,17,18]), the 
UV - soft X-ray emission from the disk — either entering the jet directly (the 
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FIGURE 1. Illustration of the model geometry and the relevant 7 radiation mechanisms for 
leptonic jet models. 

ECD [External Comptonization of Direct disk radiation] process; [13,19]) or af- 
ter reprocessing at the broad line regions or other circumnuclear material (the 
ECC [External Comptonization of radiation from Clouds] process; [14,20,21]), or 
jet synchrotron radiation reflected at the broad line regions (the RSy [Reflected 
Synchrotron] mechanism; [22-24]). 

The relative importance of these components may be estimated by comparing 
the energy densities of the respective target photon fields. Denoting by u'^ the co- 
moving energy density of the magnetic field, the energy density of the synchrotron 
radiation field, governing the luminosity of the SSC component, may be estimated 
by u'gy ~ u'b tt 7e , where tt = n'^^^ R'^ ax is the Thomson depth of the relativistic 
plasma blob and 7e is the average Lorentz factor of electrons in the blob. The SSC 
spectrum exhibits a broad hump without strong spectral break, peaking around 
(e)s5c ~ (B'/Bcr) -0 7^ ~ {^)syl1i where B' is the co-moving magnetic field, B^r = 
4.414 ■ 10^'^ G, and D = (T [1 — Pr cos9obs])~ is the Doppler factor associated 
with the bulk motion of the blob. Throughout this paper, all photon energies are 
described by the dimensionless quantity e = hu/^nieC^). 

If the blob is sufficiently far from the central engine of the AGN so that the 
accretion disk can be approximated as a point source of photons, its photon en- 
ergy density (in the co-moving frame) is m'^, ~ Ljy/lAir z"^ cT"^), where Ljj is the 
accretion disk luminosity, and z is the height of the blob above the accretion disk. 
The ECD spectrum can exhibit a strong spectral break, depending on the exis- 
tence of a low-energy cutoff in the electron distribution function, and peaks at 
{^)ecd ~ {^)d (-D/r)7g, where (e)^) is the average photon energy of the accretion 
disk radiation (typically of order 10~^ for Shakura-Sunyaev type accretion disks 
[25] around black holes of ~ 10^ - IO^^Mq). 

Part of the accretion disk and the synchrotron radiation will be reprocessed by 
circumnuclear material in the broad line region and can re-enter the jet. Since 
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FIGURE 2. Fit to the simultaneous broadband spectrum of the FSRQ 3C279 during its very 
bright 7-ray flare in September f991. See [29] for model parameters. 

this reprocessed radiation is nearly isotropic in the rest-frame of the AGN, it will 
be strongly blue-shifted into the rest-frame of the relativistically moving plasma 
blob. Thus, assuming that a fraction aBLR of the radiation is rescattered into the 
jet trajectory, we find for the energy density of rescattered accretion disk photons: 
'^'ecc ^ ttBLR r^/ (47r {r)'^BLR c); where (r) blr is the average distance of the BLR 
material from the central black hole. The ECC photon spectrum peaks around 



For the synchrotron mirror mechanism, additional constraints due to light travel 
time effects need to be taken into account in order to estimate the reflected syn- 
chrotron photon energy density (for a detailed discussion see [24]), which is well 
approximated by u'^sy ~ u'^yAT^ aBLRiRsI ^"^blr) (1 - 2T R'b/z), where ^tblr 
is a measure of the geometrical thickness of the broad line region. Similar to the 
SSC spectrum, the RSy spectrum does not show a strong spectral break. It peaks 



around {e)RSy ~ {B' / B,,) DV^ ~ {^)sscT\ 
TRENDS BETWEEN DIFFERENT BLAZAR CLASSES 



There appears to be a more or less continuous sequence in the broadband spectral 
properties of blazars, ranging from FSRQs over LBLs to HBLs, which was first 
presented in a systematic way in [26]. While in FSRQs the synchrotron and 7- 
ray peaks are typically located at infrared and MeV - GeV energies, respectively, 
they are shifted towards higher frequencies in BL Lacs, occurring at medium to 
even hard X-rays and at multi-GeV - TeV energies in some HBLs. The bolometric 
luminosity of FSRQs is — at least during 7-ray high states — strongly dominated 
by the 7-ray emission, while in HBLs the relative power outputs in synchrotron 
and 7-ray emission are comparable. 

Detailed modeling of several blazars has indicated that this sequence appears to 
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FIGURE 3. Fit to a weekly-avcragcd broadband spectrum of the extreme HBL Mrk 501 during a 
high 7-ray state centered on MJD 50564 in 1997. The dotted "bow-tie" curve indicates the highest 
flux ever measured by EGRET, which provides an upper hmit for the possible contribution of 
external Comptonization. See [32] for model parameters. 



be related to the relative contribution of the external Comptonization mechanisms 
ECD and ECC to the 7-ray spectrum. While most FSRQs are successfully modelled 
with external Comptonization models (e. g., [21,27-29]), the broadband spectra 
of HBLs are consistent with pure SSC models (e. g., [30-32]). BL Lacertae, a 
LBL, appears to be intermediate between these two extremes, requiring an external 
Comptonization component to explain the EGRET spectrum [33,34]. Figs. 2 and 
3 illustrate detailed modeling results of two objects located at opposite ends of this 
sequence of blazars, using the jet radiation transfer code described in [15,34]. 

A physical interpretation of this sequence in the framework of a unified jet model 
for blazars was given in [35]. Assume that the average energy of electrons, 7e, is 
determined by the balance of an energy-independent acceleration rate jacc and 
radiative losses, 'jrad ~ — (4/3) cctt (-u'/meC^) 7^, where the target photon density 
u' is the sum of the sources intrinsic to the jet, u'^ + u'^y plus external photon 
sources, u'^^jj^ + u'^f^Q + u'j^^y. The average electron energy will then be 7e oc 
{iacc/u'Y^"^ . If one assumes that the properties determining the acceleration rate of 
relativistic electrons do not vary significantly between different blazar subclasses, 
then an increasing energy density of the external radiation field will obviously lead 
to a stronger radiation component due to external Comptonization, but also to a 
decreasing average electron energy 7e, implying that the peak frequencies of both 
spectral components are displaced towards lower frequencies. 



SPECTRAL VARIABILITY OF BLAZARS 



Between fiaring and non-fiaring states, blazars show very distinct spectral vari- 
ability. Not only does the emission at the highest frequencies generally vary on 
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FIGURE 4. The dependence of the fit parameters F and 71 (low-energy cut-ofi' of the electron 
distribution) on the EGRET flux for fits to simultaneous broadband spectra of PKS 0528+134 
(see [28]). 

the shortest time scales, but also the flaring amplitudes are signiflcantly differ- 
ent among different wavelength bands. FSRQs often show spectral hardening of 
their 7-ray spectra during 7-ray ffares (e. g., [36,39,40]), and the flaring amplitude 
in 7-rays is generally larger than in all other wavelength bands. The concept of 
multi-component 7-ray spectra of quasars, as first suggested for PKS 0528-1-134 in 
[36] , offers a plausible explanation for this spectral variability due to the different 
beaming patterns of different radiation mechanisms, as pointed out in [37]. This 
has been applied to PKS 0528+134 in [38] and [28]. 

The results of [28] indicate that 7-ray flaring states of PKS 0528+134 are consis- 
tent with an increasing bulk Lorentz factor F of ejected jet material, while at the 
same time the low-energy cutoff 71 of the electron distribution injected into the jet 
is lowered. This is in agreement with the physical picture that due to an increasing 
r, the quasi-isotropic external photon fleld is more strongly Lorentz boosted into 
the blob rest frame, leading to stronger external Compton losses, implying a lower 
value of 7i. The external Compton 7-ray components depend much more strongly 
on the bulk Lorentz factor than the synchrotron and SSC components do. This 
leads naturally to a hardening of the 7-ray spectrum, if the SSC mechanisms plays 
an important or even dominant role in the X-ray — soft 7-ray regime, while exter- 
nal Comptonization is the dominant radiation mechanism at higher 7-ray energies. 
The results of this study on PKS 0528+134 are discussed in more detail in [41]. 

While this flaring mechanism is plausible for FSRQs, short-timcscale, correlated 
X-ray and 7-ray flares of the HBLs Mrk 421 and Mrk 501 [30,31] and synchrotron 
flares of other HBLs (e. g., PKS 2155-304, [42,43]) have been explained successfully 
in the context of SSC models where flares are related to an increase of the maximum 
electron energy, 72, and a hardening of the electron spectrum. How the spectral 
evolution in synchrotron flares can be used to constrain the magnetic fleld and 
the physics of injection and acceleration of relativistic pairs, has been described in 
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FIGURE 5. Temporal variation of the fit parameters ne(7 > 10^) (density of high-energy 
electrons) and p (spectral index of injected electron distribution) compared to the weekly averaged 
light curves from RXTE ASM, BATSE and HEGRA (see [32]). 



detail in the previous talk by R. Sambruna [44]. 

Comparing detailed spectral fits to weekly averaged broadband spectra of 
Mrk 501 [32] over a period of 6 months, we have found that TeV and hard X- 
ray high states on intermediate timcscales are consistent with a hardening of the 
electron spectrum (decreasing spectral index) and an increasing number density of 
high-energy electrons, while the value of 72 has only minor influence on the weekly 
averaged spectra. Fig. 5 shows how the spectral index of the injected electron dis- 
tribution and the density of high-energy electrons resulting from our fits are varying 
in comparison to the RXTE ASM, BATSE, and HEGRA 1.5 TeV light curves. For 
a more detailed discussion of this analysis see [45] . 

These variability studies seem to indicate that due to the different dominant 7 
radiation mechanisms in quasars and HBLs also the physics of 7-ray flares and 
extended high states is considerably different. While in FSRQs the 7-ray emission 
and its flaring behavior appears to be dominated by conditions of the external 
radiation field, this influence is unimportant in the case of HBLs where emission 
lines are very weak or absent, implying that the BLR might be very dilute, leading 
to a very weak external radiation field, which becomes negligible compared to the 
synchrotron radiation field intrinsic to the jet. 



A TOY MODEL FOR SPECTRAL VARIABILITY 

On the basis of the estimates of the photon energy densities of the various ra- 
diation fields and the peak energies of the diverse radiation components given in 
Section 2, one can develop a very simple, quasi- analytic toy model which allows us 

to study the influence of parameter variations on the predicted broadband spectrum 
of a blazar. For construction of this toy model, I assume that the magnetic field 
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FIGURE 6. Thick solid curve: Toy model calculation representative of the broadband spectrum 
of PKS 0528+134. Parameters: T = 5, dohs = 5°, 7acc = 5 • 10"* s-\ Lzj = 6 • lO'^^ erg s^S 
(e) D = 10~^ asLfi = 5 • lO^^, z^2- 10^^ cm, -R^ = 3 • IQi*^ cm, tt = lO^^. For the other curves, 
a single parameter, as indicated by the label, has been changed (value of the original calculation 
in parantheses). 



within the jet is in cquipartition with the ultrarelativistic electrons, and that the 
hght-travel time effects affecting the efficiency of the synchrotron mirror mechanism 
can be parametrized by a correction factor fiu ^ 0.1 so that v^jigy ~ v^gy (^blr^^ fitt- 
Then, the entire broadband spectrum, accounting for all synchrotron and inverse- 
Compton components, is determined by 9 parameters: F, 9obs, jacc, Ld, {^)d, clblRi 
the scale height z of energy dissipation in the jet, i?^, and tt- In most cases, several 
of these parameters can be constrained by independent observations. The radiation 
spectra of each individual component are approximated by double power-laws with 
a smooth transition. 

The thick solid curve in Fig. 6 shows a toy model calculation with the location 
of peak frequencies and the relative contributions of the 7 radiation components as 
found in our fits to PKS 0528+134 based on detailed simulations [28]. The hard 
X-ray to soft 7-ray spectrum below ~ 1 MeV is dominated by the SSC mechanism, 
while at higher energies, the ECC mechanism is dominant. The other curves in 
Fig. 6 indicate the effect of single parameter changes on the broadband spectrum 
which could be thought of as the cause of flares at 7-ray energies. 

From Fig. 6, one can see that an increasing bulk Lorentz factor leads to a strong 
flare at 7-ray energies, while only moderate flaring at infrared and X-ray frequencies 
results (dotted curve). A shift of the synchrotron peak towards lower frequencies is 
predicted [46]. If the BLR albedo oblr increases (short-dashed curve), the result is 
a slight increase in the power output at high-energy 7-rays, while due to enhanced 
external- Compton cooling the flux in the synchrotron and SSC components even 
decreases. An increased acceleration rate '^acc (long-dashed curve), leads to a strong 
synchrotron and SSC flare, where the largest variability amplitude is expected at 
MeV energies, while only moderate variability at higher frequencies is predicted. 
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FIGURE 7. Variation of the power output in synchrotron and the external Comptonization 
components if the bulk Lorcntz factor F is the only parameter changing during a flare. The 
numbers along the curves are the respective values of F at that point, and the arrows indicate 
the sense of evolution if F is increasing. 



If the density of relativistic electrons in the blob increases during a flare (dot- 
dashed curve), the variability amplitude is again predicted to be largest at MeV 
energies and the peak frequencies of all components are expected to be shifted 
slightly towards lower frequencies. A decreasing observing angle — which could be 
a consequence of a bending jet — leads to a shift of the entire broadband spectrum 
towards higher fluxes and slightly higher peak frequencies. 

Obviously, definite conclusions can not be drawn from this simplistic toy-model 
analysis. However, it may support previous results that an increasing bulk Lorcntz 
factor is a viable and plausible explanation for the spectral variability observed in 
PKS 0528+134 and possibly also in other FSRQs. 

Fig. 7 illustrates the variation of the power output in the different external 
Compton 7-ray components with respect to the synchrotron component, if F is 
the only parameter changing during a 7-ray flare. Most notably, a very steep 
relation between vFl^^^{sy) and vFl'^'^^iECC) is predicted for values of the Lorentz 
factor close to the critical F, at which the observer is looking at the super luminal 
angle. This dependence can be significantly steeper than quadratic, Ai/F^^"'^(sy) oc 

[\vFl'^°'^{ECC) with q; > 2, which has recently been observed in the prominent 

996 flare of 3C279 [40]. 
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